Regolith grain sizes of Saturn's rings inferred from 
Cassini-CIRS far-infrared spectra 



Ryuji Morishima a,b '*, Scott G. Edgington a , Linda Spilker a 

a Jet Propulsion Laboratory/California Institute of Technology, Pasadena, CA 91109, 

USA 

6 University of California, Los Angels, Institute of Geophysics and Planetary Physics, Los 
Ifi \ Angeles, CA 90095, USA 



Abstract 



We analyze far-infrared (10-650 cm -1 ) emissivity spectra of Saturn's main 
rings obtained by the Cassini Composite Infrared Spectrometer (CIRS). In 
modeling of the spectra, the single scattering albedos of regolith grains 
are calculated using the Mie theory, diffraction is removed with the delta- 
Eddington approximation, and the hemispherical emissivities of macroscopic 
free-floating ring particles are calculated using the Hapke's isotropic scat- 
tering model. Only pure crystalline water ice is considered and the size 
distribution of regolith grains is estimated. We find that good fits are ob- 
tained if the size distribution is broad ranging from 1 fim to 1-10 cm with 
a power law index of ~ 3. This means that the largest regolith grains are 
comparable to the smallest free-floating particles in size and that the power 
law indices for both free-floating particles and regolith grains are similar to 
each other. The apparent relative abundance of small grains increases with 
decreasing solar phase angle (or increasing mean temperature). This trend is 
particularly strong for the C ring and is probably caused by eclipse cooling 
in Saturn's shadow, which relatively suppresses warming up of grains larger 
than the thermal skin depth (~ 1 mm) under subsequent solar illumination. 
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1. Introduction 

Saturn's rings consist of a large number of icy particles. The range of the 
particle size in the main rings (the A, B, C rings, and the Cassini division) 
deduced from radio and stellar light occultations is roughly 1 cm to 10 m 
(Marouf et al., 1983; Zebker et al., 1985; French and Nicholson, 2000; Cuzzi 
et al., 2009). The composition of ring particles is mostly crystalline water ice 
and the mass fraction of contaminants (e.g., Tholins, PAHs, or nanohematite) 
is ten percent at most (Epstein et al., 1984; Cuzzi et al., 2009) and probably 
less than one percent (Poulet et al., 2003). A favorable origin of Saturn's 
rings with such a high content of water ice is stripping of the icy mantle of a 
Titan-sized satellite (Canup, 2010). Subsequent meteoritic pollution for 4.5 
Gyr may darken the rings too much if the ring mass has been kept similar to 
the present mass (Cuzzi and Estrada, 1998; Elliot and Esposito, 2011), but 
initially very massive rings suggested from the Canup's model have not yet 
been taken into account in pollution models. 

Individual ring particles are likely to be covered by regolith grains 0, 
probably formed by meteoritic bombardment (Elliot and Esposito, 2011). 
The regolith grain size is estimated from ring spectra at different wavelengths 
from ultraviolet to submillimeter ranges, but the estimated size is puzzling as 
it increases with wavelength. For far-ultraviolet wavelengths, Bradley et al. 
(2010) estimate the mean photon path length, which is probably the order 
of the grain size, as 2-5 fim. For near-infrared wavelengths, the estimated 
grain size varies from author to author (Nicholson et al., 2008; Cuzzi et al., 
200S|; Filacchione et al., 2012) ranging from 5 to 100 /im, but all these 
works suggest a relatively smaller size in the C ring than the A and B rings. 
These works for near-infrared spectra assume that the ring composition is 
pure water ice, and the size difference could be due to different water bands 
used or different radiative transfer models. Poulet et al. (2003) study the 
ring composition using both visible and near infrared spectra. Their model 
calls for a wide spread in grain sizes: 10-1000 \im for the A and B rings and 
30-7500 \im for the C ring. The large difference in the grain size of the C 
ring from other works indicates that there is some degeneracy between the 
ring composition and grain size. For far-infrared wavelengths, Spilker et al. 



1 Throughout the paper, we use "particles" for free-floating particles and "grains" for 
regolith grains. 

2 A work by R. Clark in Cuzzi et al. (2009) 
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(2005) find that the B ring spectrum is well fitted by a model spectrum with 
the grain sizes ranging from 8 fim to 10 m and the power-law index of 3.4 
(the upper limit of 10 m means that the entire ring is considered as regolith 
layers in their study). 

Since Saturn orbit insertion of the Cassini spacecraft in July 2004, the 
Cassini Composite Infrared Spectrometer (CIRS) has obtained millions of 
spectra (7 fim - 1 mm) of Saturn's rings (Flasar et al., 2005; Spilker et al., 
2005, 2006; Altobelli et al., 2007, 2008, 2009; Leyrat et al., 2008b; Flandes 
et al., 2010, Morishima et al., 2011). Except for Spilker et al. (2005), all 
above works discuss ring temperatures. The signal-to-noise levels of indi- 
vidual spectra are usually good enough to derive temperatures using Planck 
fits whereas averaging over many spectra are necessary for spectroscopy to 
examine grain sizes and possible contaminants. In Spilker et al. (2005), 
only the spectral data in the very early phase of the mission are analyzed. 
In the present study, using the many spectra we have obtained so far, radial 
variation and dependence on observational geometries are examined in detail. 

In Section 2, we explain how to derive emissivity spectra from observed 
radiances. In Section 3, selection and averaging of data are discussed. In 
Section 4, the modeling of emissivity spectra is described. In Section 5, the 
estimated grain size distributions for different rings are shown. Interpre- 
tations of the results and comparison with previous works are discussed in 
Section 6. The summary of the present work is given in Section 7. 

2. Derivation of emissivity spectra from observed radiances 

The quantity directly obtained by the spectrometer is the radiance 
as a function of wavenumber, v. If the emissivity e[y) is assumed not to vary 
over a footprint, the form of l[y) for ring thermal emission is given as 

I{y) = e(u)(3 gco J f(T)B(T, is)dT, (1) 

where /3 geo (< 1) is the geometric filling factor of a ring (see Sec. 2.4), f(T) is 
the distribution function of temperature T inside a footprint, normalized as 
/ f(T)dT = 1, and B(T, v) is the Planck function. Because the temperature 
inside the footprint is not known and may not be uniform, it is not straight- 
forward to derive e{y) from /(^), unlike temperature-controlled laboratory 
measurements. With different assumptions, different types of emissivities 
are obtained. Three different emissivities are introduced in the following: 
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e (v) derived by a fit with a single temperature, e\(v) derived by a fit with 
two temperatures, and £2(2^) derived by a fit with two temperatures and an 
expected emissivity spectrum. 

2.1. Fits with a single temperature 

If the temperature distribution is represented by a single effective tem- 
perature T , one may have 

J f(T)B(T, v)dT = Ahen»5(T , u), (2) 

where /3 t ^ erra is the correction factor due to multiple temperatures inside a 
footprint. The expected shape of Aherm is discussed in Sec. 2.2. 

In the first method to obtain the emissivity is to assume /3therm(^) = 1- 
Substituting Eq. (j2J) into Eq. (JTJ, we obtain 

_ I(v) 

f3 geo B(T , v) 

The temperature To is obtained by minimizing the residuals Rq 

p frp \ ST ( ^ ~ e sca,0/3gcog(T , V) V 

ito(Jo,e SC a,oJ = 2^ ^ ^/j J » I 4 ) 

where e scaj0 is the (scalar) emissivity averaged over wavelengths and a{u) 
is the instrument noise equivalent spectral radiance (NESR; see Flasar et 
al. 2004). The above Planck fits are applied to all CIRS ring spectra to 
obtain corresponding temperatures. In Spilker et al. (2005), Eqs. (J3J) and 
(J4]) are used and eo(u) is averaged over many spectra. Altobelli et al. (2008) 
estimated e scai o for the C ring to be ~ 0.9. It is expected that Eq. (j4]) works 
well if the temperature variation in the footprint is small (Ahcrm(^) — 1) and 
t{v) is roughly flat over wavelength of interest. 

2.2. Fits with two temperatures 

It is unlikely that the temperature inside a single footprint is perfectly 
uniform. The particle temperature may vary radially and vertically. Even 
on the surface of a single particle, large temperature variation may exist as 
the temperature is expected to be highest and lowest around the sub-solar 
and anti-sub-solar points, respectively. Instead of deriving a temperature 
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distribution with complicated thermal models which calculate the energy 
balance of rings, simply two portions with high and low temperatures (T wl 
and T c i) are considered: 

f(T)B(T, v)dT = f wl B(T wl , v) + (1 - f wl )B(T cl , u), (5) 

where f w i is the fraction of the warm portion. Fits with two temperatures 
are also used in analysis of thermal emission of icy satellites (Carvano et al. 
2007; Howett et al. 2011). Substituting Eq. © into Eq. flD, the emissivity 
is given by 

ei(Z/) = /3 g eo [UiB(T wU + Ui)B(T cU v)} ■ (6) 

As in Eq. 01]), the temperatures and factors are determined by minimizing 
the residual as 



Ri[T w i, T c i, e sca) i, 



wl; 



E 



~ e S ca,l/3geo [/wl-B(T w i, v) + (l - f wl )B(T cl , u)] 



<T[V) 

V 
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where e sca) i is the scalar emissivity. For the convenience of later use, we also 
define the relative emissivity e rc i^(u) as 

ereuM = (8) 
^sca,l 

The relative emissivity is roughly scaled to be unity, so is suitable for com- 
parison between different spectra. 

Equating Eqs. (j2J) and (jSJ), the form of /3therm is given as 

a f \ UiB(T wl ,u) + (l-f wl )B(T cl ,u) 

Ptherm^J = s • (9) 

An example of Ahcrm for various values of / w i is shown in Fig. 1. Near 
the Planck peak (~ 200 cm -1 ), /3 t herm takes the lowest value whereas Ahcrm 
is much larger than unity at large v. Indeed most of spectra show large e 
at large v due to the effect of /3therm, if fits are done with a single Planck 
function. 
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2.3. Fits with two temperatures and an expected emissivity 

The method in Sec. 2.2 works only if e(u) is globally flat with possible 
narrow features in the wavenumber range of interest. If broad features exist, 
correct temperatures may not be retrieved. An extreme case is that e(p) 
has a profile similar to l//3therm(^)- In this case, both e(v) and /3therm(^) 
are incorrectly found to be nearly constant over the range of v. In fact, we 
have this problem because eiv) for water ice is an upward convex function 
around the Planck peak wavenumber for the ring temperature (Sec. 4). In 
this kind of situation, the temperatures (T w2 and T c2 ) and the fraction of the 
warm portion (/ W 2) may be retrieved with an expected (modeled) emissivity 
e m odci(^i, n2, v), where n%, n 2 , ... are parameters in models, as 



R2(T w 2i T c2 , / W 2, ni, n 2 , ...) — ^ ( 

1/ ^ 



2 

— lm_ Ac\(Tw2-, T C 2, fw2, "-1,^2, V , 



,(io) 



where 



^model(^w2j 7c2 ; /w2 ; n l, n 2i •••> v ) — 

e m od e i(ni,n 2 , z/)/3 geo [/ W 25(T w2 , v) + (1 - f w2 )B(T c2 , v)} . (11) 
The retrieved emissivity is given as 

e2(Z/) = /3 geo [/w 2 5(T w2 , v) + (1 - f w2 )B(T c2l u)] ' (12) 

It is expected that e m0( j e i(^) with best-fit parameters is close to e 2 (z/). 

In this study, we estimate best-fit values of model parameters using 
Eq. (jTUl) . rather than comparing e\{v) with e mo ^ e \{v) : because only poor fits 
are obtained for the latter case. 

2-4- Geometric filling factor 

The geometric filling factor /3 geo is given as 

/W£o) = l-exp(-^|V), (13) 
x sin a 1 ' 



o I 



where B Q is the elevation angle of the observer and t(B d ) is the normal 
optical depth of the ring. The optical depth r is known to depend on B 
and also on longitude of the observer around the footprint relative to the 
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ring longitude because of wakes (Colwell et al., 2006, 2007, 2010; Hedman 
et al., 2007). In the present study, we only consider the _B -dependence and 
t(B ) is obtained by interpolations using five r profiles with different values 
of B Q (15.3, 21.1, 28.7, 54.0, and 66.7 degs). These data are available in 
Planetary Data System (Colwell et al., 2010). The r profile for B Q = 28.7 
degs is obtained by Voyager PPS while others are by Cassini UVIS. We 
interpolate r if t(B = 28.7 degs) < 0.2 then convert it to /3 geo , whereas 
/3 geo is interpolated for larger r. Usually r has little dependence on B D in 
regions with small r. The interpolated profile of /3 geo is smoothed over CIRS 
footprints (see Morishima et al., 2010). 

3. Selection and averaging of data 

Individual CIRS spectra are usually too noisy for spectroscopy and aver- 
aging over many spectra is necessary to reduce the noise level. Different types 
of averaging are considered. In previous works (Spilker et al., 2005; Carvano 
et al., 2007), the emissivities (e (^) and ei(z/)) are averaged over after ap- 
plying Planck fits to individual spectra. On the other hand, the radiance 
is averaged over and Planck fits are applied to the averaged spectra in this 
study. There are advantages and disadvantages in both methods. The latter 
method is numerically less intense than the former method, particularly if 
the number of parameters used in Planck fits are numerous like the case of 
e 2 (z^). The channel-to-channel noise level is smaller for the latter method 
than the former method, as the noise originally arises in the radiance. 

However, the following conditions are necessary for the latter method to 
work: (1) variation of the emissivity is small enough in the spectra used in 
an average, (2) correlation between /3 therm _B(T, u) (temperature and its vari- 
ation) and /9geo(-B ) is small enough, and (3) temperature variation including 
both inside of individual footprints and over different footprints are still well 
approximated by a sum of two Planck functions. Averaging Eq. ([1]) over 
spectra, these conditions are represented as: 

(I(y)) = e 1 (z/)(Ahcrm(^)5(T,z/)/3 geo ) 
= ei(v)(f3 thcTm (v)B(T,v))((3 geo ) 

= £!(» [f wl B(T wl , u) + (1 - Ux)B(T cl: u)\ (/3 geo ). (14) 

A similar equation is also obtained for ^{v)- 

In order to check whether these conditions are fulfilled, we test various 
types of grouping of spectra based on saturnocentric radius, observational 
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geometry, and temperature, T , while keeping a certain number of spectra 
in a parameter bin for a good signal-to-noise level. After many trials, we 
divide spectra into six radial bins (the inner and outer C ring, inner and 
outer B ring, Cassini division, and A ring) and into temperature (T ) bins 
with the interval of 5 K. These bins are chosen because apparent spectral 
shapes depend on saturnocentric radius and temperature. The temperature 
dependence probably represents the dependence on observational geometries. 
Noise cancellation works better with temperature bins rather than using bins 
of geometric parameters, and also an unnecessary enhancement of the effect 
of Ahcrm is suppressed by averaging spectra with similar temperatures. We 
do not divide spectra by B Q because we find that there is only small difference 
between averaged spectra at high and low \B Q \, provided that other condi- 
tions are the same. This means that the condition (2) is almost fulfilled. 
The number of spectra in the bins, the mean solar phase angle, and the 
mean absolute solar elevation angles are shown in Fig. 2. The highest tem- 
perature is obtained at low phase angle and at high solar elevation, while the 
temperature decreases with decreasing solar elevation and increasing phase 
angles. 

The above binning is made for all spectra obtained until the end of 2010 
after removing very noisy spectra. We only analyze spectra in the far-infrared 
channel (focal plane 1; 10 cur 1 < v < 650 cm -1 ) with T > 70 K, as spectra 
at low temperatures or in the mid-infrared channel (focal planes 3 and 4; 600 
cur 1 < v < 1400cm -1 ) are too noisy even after averaging. We only show 
spectral data at the lowest spectral resolution, 15.5 cm -1 (the interpolated 
spectral step is 5 cur 1 ), in this paper. Spectra with higher resolutions are 
noisier but show good agreement with the lowest resolution data in most 
cases. 

Examples of the relative emissivity e rc i ) i(z/) for different temperatures T 
are shown in Figs. 3 and 4 for the inner C ring and the outer B ring, respec- 
tively. The emissivity spectra for the inner C ring show strong dependence 
on temperature whereas the inner B ring spectra show little dependence on 
temperature. These temperature dependence and independence are also con- 
firmed in higher spectral resolution data. The dip seen at 230 cm -1 (clearly 
seen in Fig. 4 but ambiguous in Fig. 3) is a feature caused by intermolecular 
vibrations (transverse-optical mode) of water ice. Other narrow features are 
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noises specific to the CIRS instrument (Carlson et al. 2009j): the real-time 
interrupt (RTI) noise (191, 383, and 574 cm" 1 ) and a prominent single fre- 
quency (sine wave in interferograms) feature of unknown origin that has been 
observed between 18 cm -1 and 210 cm -1 (seen at 80 cm -1 in Fig. 3 and at 
150 cm -1 in Fig. 4). The noise level is not shown in Figs. 3 and 4, but the 

1 /2 

channel-to-channel noise level is expected to be ~ a(y) /((f3 gco )B(T , ^)iV S pe C ) 
and is small enough due to large N spec , except at the wave numbers of the 
CIRS specific noise. The standard deviation of the radiance in spectral aver- 
aging is not used, as it is independent of iV S p ec and makes formal error bars 
artificially large (Clark et al., 2008). Most of the spectra in Figs. 3 and 4 
show increase of e re i,i(z/) with v for v > 400 cm" 1 . This is due to the effect 
of Aherm which is not removed appropriately even with two temperature fits. 
This artificial effect will be better corrected in ^iy)- 



4. Emissivity modeling 

In our emissivity model, the single scattering albedos of regolith grains 
are calculated using the Mie theory (Mie, 1908; Bohren and Huffman, 1983), 
diffraction is removed with the delta-Eddington approximation (Joseph et al., 
1974), and the hemispherical emissivities for macroscopic particles are cal- 
culated using the Hapke's isotropic scattering model (Hapke, 1993). Cheng 
et al. (2010) compare the mid-infrared emissivity spectra (8-13 /xm) of snow 
surfaces obtained by field measurements of Hori et al. (2006) with those from 
various types of emissivity models, and find that the Mie-Hapke hybrid model 
with diffraction removal is the best choice. The Mie-Hapke hybrid model is 
also the best choice for emissivity modeling of quartz (Moersch and Chris- 
tensen, 1995) and olivine (Mustard and Hays, 1997), although diffraction 
removal does not improve fits for quartz (Pitman et al., 2005). 

The Mie theory calculates absorption and scattering of a sphere with a 
given size parameter X = 2nr/X (where r is the grain radius and A is the 
wavelength of incident wave) and the real and imaginary parts of the refrac- 
tive indices, n and k. In the present work, we only consider pure crystalline 
water ice. To date non-water ices (NH3, CO2, and CH4) have not yet been 



3 De-spiking of interferograms described in Carlson et al. (2009) has not been applied to 
the 15.5 cm -1 data. See also "Interferences on CIRS interferograms and spectra: A user 



guide" by Nixon et al. at |http://pds-rings.seti.org7 ; TO l/COCIRS-0409/DOCUMENT/cirs. 
intcrfcrences.pdf. 
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identified in CIRS far-infrared spectra (Edgington et al., 2008), and effects of 
possible contaminants will be examined more in detail in future work. The 
refractive indices of crystalline water ice used for the Mie calculations are 
shown in Fig. 5. We use n and k at 136 K from Curtis et al. (2005) for 50 
cm -1 < v < 650cm -1 and k at 90 K from the theoretical fitting function 
of Mishima et al. (1983) for 10 cm" 1 < v < 30 cm . We assume that n 
is constant for v < 50 cm -1 , and log A; is linearly interpolated between 30 
and 50 cm -1 . The temperature of Saturn's rings is lower than 136 K, but no 
good data are available at ring temperatures for crystalline water ice. The 
refractive indices weakly depend on temperature, particularly, v at the peak 
of k (near 230 cm -1 ) weakly increases with decreasing temperature (Curtis 
et al., 2005). In fact, the CIRS emissivity spectra show the peak wavenumber 
slightly larger than the model prediction (see Fig. 8), as we use n and k at the 
warmer temperature. The samples of Curtis et al. (2005) at temperatures 
lower than 136 K are amorphous and amorphous ice shows u at the peak 
of k slightly lower than that for crystalline ice (Fig. 5). This indicates that 
Saturn's rings most likely consist of crystalline ice. The same conclusion is 
also derived from the Fresnel peak position in near-infrared spectra (Cuzzi 
et al., 2009). 

Using the extinction efficiency Q cx to and the scattering efficiency Qscao 
from the Mie theory, the single scattering albedo u is given as 

QscaO ; 1r \ 

= 75 — • (15) 

VcxtO 

Figure 6a shows Q cx to and u for various grain radii. In the Mie theory, 
scattered light includes diffracted light. In a closely packed medium, however, 
inter-grain spaces are narrow so that the diffracted light is suppressed (Hapke, 
1999). Different methods of diffraction removal are proposed (Joseph et al., 
1976; Wiscombe, 1977; Mishchenko, 1994; Wald, 1994) and we apply the 
delta-Eddington approximation proposed by Joseph et al. (1976) because 
of its convenience. In this approximation, the diffraction is approximated 
by a forward-scattering delta function and its contribution is calculated by 
assuming that the second moment of the sum of the delta function and the 
diffraction-removed phase function is identical to that for the original phase 
function. The extinction efficiency Q e xt, the scattering efficiency Q sca , and 
the single scattering albedo, u, after diffraction removal are given as 

Qext = (1 - UogDQextO, (16) 
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Q, 



sea 



(1 ~~ 9o)^oQscaX), 



(17) 
(18) 



sea 



CO 



Qcxt 



where g$ is the anisotropic parameter calculated by the Mie-theory. Fig- 
ure 6b shows Q ext and u for same grain sizes used in Fig. 6a. It is found that 
oj is usually much lower than ujq. For X 1, Q cx to — 2 whereas Q cxt — 1, 
as expected. The Mie parameters are little affected by the delta-Eddington 
approximation for small grains (X < 1), because go approaches zero in the 
Rayleigh scattering regime. However, differences in u due to different diffrac- 
tion subtraction methods are large for small grains (not shown). Fortunately, 
this has little impact if a broad size distribution including large grains is taken 
into account, because Qext for small grains is so small that its contribution 
to scattering is reduced (see Eq. (fl9|) ). The exception is near the peak of k 
at 230 cm -1 , where Q ext is large even for X < 1, and there is little difference 
between Q ex to an d Qext, even applying any diffraction subtraction method. 
For a size distribution, the single scattering albedo is effectively given as 



We assume the number of grains per unit radius is given by a power-law 



where p is the index, n and r are constants. The black thin lines in Fig. 6a 
and 6b show Uq and u for p = 3 with the minimum and maximum cut-off sizes 
of r m ; n = 1 /im and r max = 1 cm. The logarithmic bins with an increment of 



Hapke (1993) derives analytic expressions of various types of reflectances 
and emittances using the two stream approximation for isotropic scattering 
grains. Using u, the Bond (spherical) albedo of a macroscopic ring particle 
comprised of regolith grains is given as (Eq. (10.51b) of Hapke, 1993) 



u = 



J 7ir 2 Q sca (r)n(r)dr 
J 7rr 2 Q ext (r)n(r)dr' 



(19) 




(20) 



2 1//4 are used for r. 





with 



7 = vi — UJ. 



(22) 
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Finally, the hemispherical emissivity e mo dc\ of the ring particle is given from 
the Kirchhoff's law as (Eq. (13.29) of Hapke, 1993) 



For spatially resolved planetary or satellite surfaces, the directional emissiv- 
ity is usually used, whereas we use the hemispherical emissivity (the hemi- 
spherically averaged directional emissivity) as we see thermal emission from 
entire hemispheres of ring particles. We only considered intra-particle re- 
golith radiative transfer, ignoring inter-particle scattering. Indirect evidence 
that inter-particle scattering is unimportant is obtained from a similarity of 
the emissivity spectra between the B and A rings despite the optical depth 
difference between them (see Sec. 5 for a more detailed discussion). 

Bottom panels of Figs. 6a and 6b show e mo dei without and with diffraction 
removal. For large grains (r > 100 /im), the overall profile of e mo dei is flat 
down to ~ 50 cm -1 , and a spectral roll-off is seen for lower v\ the exact 
value of v below which a roll-off is seen decreases with increasing r. The 
location of a roll-off seen in CIRS spectra (Figs. 3 and 4) indicates that large 
grains indeed exist in regoliths of Saturn's rings. The prominent dip of e mo dei 
corresponding to the k peak is seen at 230 cm -1 , and e mo dci is slightly larger 
for > 230 cm -1 than for < 230 cm -1 (for example, compare values at 170 
cm -1 and 280 cm -1 ). Intermediate-sized grains (10 /im < r < 100 /im) 
show a broad hump between 100 cm -1 and 300 cm -1 due to inter-molecular 
vibrations of water ice. The hump features are also seen in many CIRS 
spectra. The emissivity of small grains (r < 10 /im) are almost unity for 
< 300 cm -1 , because they behave as Rayleigh absorbers. Their extinction 
efficiencies are high only near 230 cm" 1 so that they reduce the depth of the 
dip at 230 cm -1 produced by large grains. 

The one issue we find in emissivity fits is that a large amount of small 
grains is necessary to produce the shallow 230 cm -1 feature seen in CIRS 
spectra whereas too large a fraction of small grains makes entire fits worse. 
An alternative way to reduce the depth of the 230 cm -1 feature instead of 
adding too many small grains is to introduce intra-grain pores. If the pore 
size is sufficiently smaller than the wavelength, the complex permittivity, 
K{4>) = (n 2 — k 2 ) + i(2nk), as a function of the porosity is given from the 
Maxwell-Garnett theory (Maxwell-Garnett, 1904) as 



Cmodel — 1 ~ ^BOND- 



(23) 



K{4>) = K(0) 1 + 



30(1 -K(0)) 



) 



(24) 



1 + 2K(0) - 0(1 - K(0)) 
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The validity of this equation is confirmed by laboratory measurements of 
silicon (Labbe-Lavigne et al., 1998). The refractive indices for = 0.5 are 
shown in Fig. 5. Figure 6c shows Q ex t, oj, and e mo d e i for (f> = 0.5 and with 
diffraction removal. It is found that the overall shape of e mo dei in Fig. 6c 
is similar to that in Fig. 6b except the depth of the 230 cm -1 feature is 
shallower with intra-grain pores. Since introducing intra-grain pores has 
an effect similar to reducing grain-size, the estimated grain size effectively 
increases, but the change in estimated size is very small as long as < 0.5. 

5. Estimated regolith size 

The best-fit values of model parameters are estimated from Eq. fllOp . In 
the model fits, we fix the minimum grain size and intra-grain porosity at 
^min = 1 /^rn and <fi = 0.5 and seek the best-fit values of the maximum 
size r max and the power-law index p. The reason for choosing r m j n = 1 /zm 
is that existence of micron-sized grains is supported from the grain size of 
spokes (r ~ 2 /im; D'aversa et al., 2010) and from UVIS spectra (r ~ 2-5 
/im; Bradley et al., 2010). We do not examine cases with r max larger than 
1 m because this is a typical size of large free-floating particles (French and 
Nicholson, 2000). 

The residual R 2 is calculated as a function of r max and p, each pair of 
r max and p having the optimal temperatures and the fraction of the warm 
portion. We introduce the normalized residual: 

i?2(p,7w) = -^--l, (25) 

-'*'2,min 

where -R2,min is the minimum value of R 2 in the parameter space. Some 
examples of R 2 on the r max vs. p plane are shown in Fig. 7. The error bars 
for r max and p are estimated from the maximum and minimum values of these 
parameters on the curve of R 2 = 1. It is found that p can be constrained 
in some ranges whereas only the lower limit of r max (3-10 mm for cases in 
Fig. 7) is well constrained. 

Some examples of the best-fit spectrum, e mo d e i, are shown together with 
e 2 and e\ in Fig. 8. Note that e 2 and e\ are obtained from the same radiance 
data and the difference comes from the fitted temperatures and the fraction 
of the warm portion (Sections 2.2 and 2.3). Good agreements between e mo doi 
and e 2 are obtained for all cases at least for v < 300 cm -1 and sometimes 
over all the wavenumber range. On the other hand, the absolute value of e x 
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is much lower than e mo d c i with any grain sizes. In addition to that, t\ for v ~ 
170 cm -1 is comparable to or higher than e\ for v ~ 280 cm -1 whereas the 
opposite trend is seen for e mo dci- The temperature T c2 is found to be always 
lower than T c \. With decreasing temperature from T c \ to T c2 , the absolute 
value of e increases in all wavenumbers and the overall slope dejdu increases. 
Consequently e 2 for v ~ 280 cm -1 can be larger than that for v ~ 170 cm -1 
as is the case of pure water ice. Although lowering only the temperature 
of the cold portion makes e at large v (typically > 400 cm -1 ) larger than 
unity e in this region is reduced by increasing the temperature and fraction 
of the warm portion (note / w2 > / w i or /and T w2 > T w \ in all cases). Some 
cases shows only one temperature component (/ w i = 0) for the derivation 
of ei (e.g., the C ring data in Fig. 8), but existence of two temperature 
components, as is the case for e 2 , is also physically reasonable considering a 
possible temperature variation over the surfaces of ring particles. 

Figure 9 shows the summary of estimated values of p and r max for all 
cases. It is found that p ~ 3 for most of the cases (Fig. 9a) except that the 
C ring shows slightly larger values and the Cassini division shows slightly 
smaller values at high temperatures. Interestingly, the values of p are quite 
close to those for macroscopic particles (French and Nicholson, 2000). The 
maximum size r max is estimated to be 1-100 cm in most of cases. The lower 
limit of r max is well constrained, and in almost all cases r max is larger than 
1 mm (Fig. 9b). On the other hand, the upper limit of r max is not well 
constrained from our data. However, r max is expected to be smaller than the 
size of large particles ~ 1 mat least by an order of magnitude and probably 
even smaller in the following reasons. First, the total cross sections of any 
particles are equal in logarithmic bins for p = 3, and thus the broad size 
distribution we obtained means that the surface of a free-floating particle is 
not occupied by the largest regolith grains. Second, if inter-grain scattering is 
less than that expected from the radiative transfer model due to large grains 
comparable to the particle in size, the emissivity spectra may become closer 
to 1 — uj instead of 1 — Aeond- This clearly makes our fits worse; in particular, 
the depth of the 230 cm -1 feature will be too deep. Most likely, r max is ~ 1 
-10 cm. This means that the largest regolith grains are comparable to the 
smallest free-floating particles in size. It is likely that grains smaller than 
1-10 cm adhesively stick to large particles due to the surface energy (Bodrova 
et al., 2012), while larger grains are cores of particles hidden under smaller 
grains, if they exists. 

The A and B rings show quite similar values of p and r max . In addition 
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to that, these rings show very similar depths of the 230 cm -1 feature. If the 
optical depth of a ring is large enough and inter-particle scattering occurs, the 
emissivity of the ring may be calculated from Eqs. ( 12 ip and (T23]) assuming 
that the single scattering albedo is now given as the Bond albedo of ring 
particles. The resulting ring emissivity should have a shallower depth of the 
230 cm' 1 feature as compared with that without inter-particle scattering, 
from comparison of the depths for 1 — u and e mo d e i in Fig. 6b or 6c. Therefore, 
if regolith sizes are the same and inter-particle scattering is important, the 
depth of the 230 cm -1 feature is expected to decrease with increasing optical 
depth; such a correlation is not seen in spectra (r ~ 2 for the outer B 
ring while r ~ 0.5 for the A ring). Therefore, inter-particle scattering of 
far-infrared light is likely to be negligible. Little inter-particle scattering in 
near-infrared light is also suggested at low phase angles (Cuzzi et al., 2009). 

The observed spectra for the C ring at low temperatures or the Cassini 
division do not clearly show the 230 cm -1 feature in contrast to the modeled 
ones. Higher spectral resolution data (not shown) seem to show that the 
feature actually exists both for the C ring and the Cassini division but it 
is narrower and shallower than that for the A and B rings (see also Figs. 3 
and 4). What causes the feature less evident is unclear. It may be due to 
low signal-to-noise ratios for the optically thin rings, to intra-grain porosities 
higher than we assume (0.5), or to effects of contaminants including amor- 
phous water ice. Since higher spectral resolution data are noisy, more careful 
analysis is necessary. 

6. Discussion 

6.1. Dependence of size distribution on temperature 

What causes the dependence of size distribution on temperature seen 
in Fig. 9? This is probably caused by effects of observational geometries 
in conjunction with non-uniform temperature distribution of ring particles 
over surfaces and depth. If the temperature of a ring particle is completely 
uniform, it is expected that the emissivity spectrum of the particle is inde- 
pendent of observational geometry. However, in actual rings, the strongest 
thermal emission comes from near the subsolar points (seen at low solar phase 
angles) of slowly rotating particles at high solar elevation angles. In addition 
to that, due to eclipse cooling in Saturn's shadow, the strongest emission 
comes from only a thin surface layer with a thickness given by the thermal 
skin depth £. If the thermal inertia is ~10 Jm -2 K -1 s -1 / 2 as estimated by 
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Ferrari et al. (2005) and Morishima et al. (2011), £ ~ 1 mm (e.g., Morishima 
et al., 2009). This is smaller than the size of the largest regolith grains es- 
timated in Sec. 5. Therefore, the largest grains are only partially heated or 
their temperatures are lower than small grains near the subsolar point. As a 
result, thermal emission from the largest grains is suppressed relative to that 
from small grains at high temperatures. At higher solar phase angles and at 
lower solar elevation angles, variation of the particle temperature with depth 
is much smaller so that large grains can have temperatures similar to those 
for small grains. The dependence of size distribution on temperature is seen 
most strongly in the C ring, probably because the periodic temperature vari- 
ation due to eclipse cooling is particularly large for the C ring as it is close to 
Saturn (long duration in the shadow) and has low photometric albedo (high 
temperature outside the shadow). 

To reinforce the above discussion, let us estimate the depth from which 
thermal emission comes. The electric skin depth may be simply given by a 
depth where the optical depth of regolith layers is unity, unless grains are 
highly reflecting. The regolith optical depth as a function of depth is given 
as 

/'''max 

t(z) = z irr 2 Q cxt (r,is)n(r)dr, (26) 



mm 



where n is the particle number density per unit volume and per unit size. 
The number density n is defined by the following equation: 



1 - 



''inter 



/ -nr 3 n(r)dr, (27) 

J r *J 



where Winter is the inter-grain porosity. Using these equations, the electric 
skin depth is given as 

z{T = V = Jn — n cg 2 v ^ 

4 Voxt,eflfU — Winter J 

where the effective grain radius r cS and the effective extinction efficiency 
Qcxt.cff are given by 

J r . nr 2 n(r)dr 

and 

/ r r ^ x 7rr 2 Q cxt (r, v)n(r)dr 

Q ex t,es{v) = mm f r max 2 , ■ (30) 
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Since Q e xt,es( l> ) is ~ 0.5-1.0, z{t = 1) is larger than r eff by a factor of 2-3, 
unless 0mter is close to unity. Figure 10a shows r e g- estimated from p and r max 
from Fig. 9. It is found that r c s is ~ 1 mm (~ £) at the highest temperatures 
and is one or two orders of magnitude larger for lower temperatures. This 
indicates that if the temperature is uniform as is likely to be the case at 
low temperatures, thermal emission come from 1-10 cm in depth. However, 
thermal emission at the highest temperatures comes from a surface layer of 
1 mm because only the thin layer is heated to the the highest temperature 
and emission from deeper and colder layers contribute much less to the total 
emission from the surface. 

6.2. Comparison with grain sizes estimated in previous works 

It is not straightforward to compare the size distribution estimated from 
the present work and those from other works, because other works usually 
consider single sizes or a few discrete sizes. If spectral behavior is monoton- 
ically dependent on grain size, a spectrum for a broad size distribution may 
be similar to the one for the median grain radius r mc d, which is defined so 
that the total cross section of particles with r > r med is equal to that with 
r < r med : 



Figure 10b shows r mcd . For the C ring, r med strongly depends on temperature, 
but as we discussed in Sec. 6.1, the relative fraction of small grains are likely 
to be overestimated at high temperatures. Therefore, r med is typically a few 
hundred microns or larger for all rings. 

The mean photon path length (2-5 fim) estimated from far-ultraviolet 
spectra (Bradley et al., 2010) is clearly much smaller than r med . The path 
length probably represents the size of cracks or flaws inside large grains or 
represents the smallest grain size, as suggested by Bradley et al. (2010). 
The typical value of r mcd seems to be also larger than the grain sizes (< 
100 jim) estimated from the water band depths in near-infrared wavelengths 
(Nicholson et al, 2008; Cuzzi et al., 2009; Filacchione et al., 2012). It is 
cautioned that this comparison is not very precise, as some of the water 
band depths do not depend on particle size monotonically (see Fig. 15.27 of 
Cuzzi et al. (2009)). Nevertheless, some other bands show rather monotonic 
dependences on size (up to ~ 1-10 cm) and the estimated sizes from these 
bands are still smaller than 100 \im. This size is probably underestimated 
as discussed in the following. 




(31) 
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The grain sizes of snow on Earth are also estimated from reflectances in 
different bands in visible and near-infrared wavelengths and the estimated 
sizes are compared with the sizes from in-situ measurements. Some works 
(Fily et al, 1997; Aoki et al, 2007; Kuchiki et al, 2009) find that the grain 
sizes estimated from the 1.6 /zm band, which is in one of the water absorption 
bands, are underestimated. Kuchiki et al. (2009) find that the grain size is 
underestimated when the actual grain size is large (> 1 mm) and point out 
that it is probably caused by microstructures on large grains' surfaces. It is 
plausible that regolith grains in ring particles are also not spheres but have 
surface microstructures. If it is the case, it may not be very appropriate 
to estimate the sizes of large grains from the water band depths, whereas 
spectral roll-offs at long wavelengths seen in CIRS spectra are likely to be 
less affected by microstructures. The large C ring grain sizes (30-7500 fim, 
which gives r mc d ~ 1 mm) estimated from spectra in the entire visible and 
near-infrared wavelength range (Poulet et al., 2003) look consistent with our 
results. 

7. Summary 

In this study, the far-infrared emissivity spectra of Saturn's main rings 
obtained by Cassini CIRS are analyzed. Ring spectra are divided into six 
radial bins and temperature bins with the width of 5 K, and the radiances 
are averaged over bins to reduce the noise levels. The temperatures and size 
distribution of regolith grains are simultaneously retrieved by the chi-square 
fits, provided that ring particles are composed of pure water ice. In modeling, 
the Mie-Hapke hybrid model is used with diffraction removal. Good fits are 
obtained if size distribution is broad ranging from 1 fim to 1-10 cm with 
the power law index of ~ 3. This means that the largest regolith grains are 
comparable to the smallest free-floating particles in size and that the power 
law indices for both free-floating particles and regolith grains are similar to 
each other. 

The apparent relative abundance of small grains increases with decreasing 
phase angle (or increasing mean temperature). This trend is particularly 
strong for the C ring and is probably caused by eclipse cooling in Saturn's 
shadow, which relatively suppresses warming up of grains larger than the 
thermal skin depth (~ 1 mm) under subsequent solar illumination. The A 
and B ring spectra are found to be very similar to each other despite the 
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difference in optical depth, and this probably indicates that there is little 
scattering of far-infrared light between free-floating macroscopic particles. 

There are a lot of things to do in future work. Effects of possible con- 
taminants need to be investigated. Since the imaginary part of the refractive 
indices for water ice is high in the broad range of far-infrared wavelengths, 
effects of contaminants are probably small if they are mixed as small in- 
clusions (r <C A) inside grains. Nevertheless, some possible contaminants, 
such as SiC>2, Fe2C>3, NH3, and NaCl, have strong absorption peaks in the 
CIRS wavelength region, and the upper limits of their abundances may be 
well constrained even if they are not identified. The effects of contaminants 
are probably more important if they are intimately mixed (r > A). One of 
plausibly intimately-mixed contaminants is amorphous carbon (Poulet et al., 
2003). It is also necessary to reduce CIRS specific noises with better cali- 
bration. Use of high spectral resolution data may also help improve model 
fits. The accurate determination of refractive indices of crystalline water ice 
at low ring temperatures (< 100 K) are also needed. 
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Figure captions 

Fig. 1. Examples of Aherm as a function of is, given by Eq. (Q. Here we 
adopt T wl = 105 K and T cl = 80 K, and T is obtained from a Planck 
fit (Eq. (j3])). Different colors represent different values of the fraction 
of the warm portion, / wl . 

Fig. 2. Number of spectra, N spec , mean solar phase angle, a, and mean 
absolute solar elevation angle \B'\ of spectrum divided in six radial 
bins and temperature (T ) bins of a 5 K interval. 

Fig. 3. The relative emissivity e re i i for the inner C ring given by Eq. (jSJ). 
Different colors represent different temperature bins: 70-75 K(purple), 
75-80 K(blue), 80-85 (light blue), 85-90 K (brown), 90-95 K (orange), and 
95-100 K (red). The upper panel shows the entire wavenumber range 
of the CIRS far-infrared channel, whereas the bottom panel shows the 
short wavenumber roll-off region. Note the difference in the vertical 
scales in these panels. 

Fig. 4. The same with Fig. 3, but for the case of the outer B ring. 

Fig. 5. The real (n) and imaginary (k) parts of refractive indices of water 
ice. The black solid lines represent n and k for crystalline water ice at 
136 K for v > 50 cm - 1 (Curtis et al, 2005) and k at 90 K for 10 cm" 1 
< v < 30 cm ~ 1 (Mishima, 1983). n is assumed to be constant for v < 
50 cm _1 , and log/c is linearly interpolated between 30 cm -1 and 50 
cm -1 . The black dashed lines represent n and k with the intra-grain 
porosity of 0.5. The blue lines represent n and k for amorphous water 
ice at 106 K (Curtis et al, 2005). 

Fig. 6. Extinction efficiency, Qext, single scattering albedo, u, and emissiv- 
ity, e mo dei, as a function of wavenumber calculated from the Mie theory 
and the Hapke's model. Different colors represent different grain sizes. 
In the panels of u and e mo d e i, the black lines represent a size distribution 
case with r min = 1 //m, r max = 1 cm, and p = 3. (a) parameters calcu- 
lated without diffraction removal (Q e xt0) and e mo dei calculated with 
uq), (b) parameters calculated with diffraction removal, and (c) param- 
eters calculated with diffraction removal and the intra-grain porosity </> 
of 0.5. For Figs. 6a and 6b, = 0. 
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Fig. 7. Contour of R 2 on the p vs. r max plane. A cross represents the 
location of the best-fit values where R 2 = 0. 

Fig. 8. Emissivities from observations e\ (blue) and t 2 (black) and modeled 
emissivity e mo d e i (red) with the best fitted r max and p (shown in the 
parenthesis) for the inner C ring, the outer B ring, the Cassini division, 
and the A ring. For each ring, one high temperature case and one low 
temperature case are shown. The fitted temperatures and the fraction 
of the warm portion are shown in the panels. 

Fig. 9. Estimated p and r max for all temperatures and all rings. 

Fig. 10. Effective radius r cS (Eq. (I2"9"j) ) and median radius r med (Eq. ( 1HT|) ) 
calculated from p and r max in Fig. 9. 
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(a) With diffraction 
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(c) Diffraction removed, = 0.5 
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nner C ring, 75-80 K 





3.4 




















______ 10.0 — > 


:' n d e x , 


3.2 








___ 5.0 = 

— 3.0 — 


o 

T 

<D 


3.0 











O 


2.8 










O. 














2.6 






— 1.0 





10 3 10 4 10 5 10 6 

Maximum size ^ttl) 



nner C ring, 95-100 K 



3.4 







^ .3.U- 1 — ~ ' 




1 .0 — - 






/ " ><_, 










- .3 — 


^ 3.0— ' 

_ 5.0 — 




7' 

v 


/ 






__io.o— : 



10 3 10 4 10 5 10 6 

Maximum size (fini) 



Outer B ring, 90-95 K 




10 3 10 4 10 5 

Maximum size (/xm) 



Fig. 7. 



36 



Inner C ring, 75 — 80 K 
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Outer B ring, 75-80 K 
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